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Primordial black holes, with masses comparable to asteroids, are an attractive possibility for 
dark matter. In addition, other forms of dark matter could form compact dark objects (CDO). We 
search for small tidal accelerations from low mass black holes or CDOs orbiting near the Earth, and 
find none. Using about 10 years of data from the superconducting gravimeters in the Black Forest 
Observatory in South-Western Germany and at Djougou, Northern Benin in Western Africa we set 
an upper limit on the maximum mass of any dark object orbiting the Earth as a function of orbital 
radius. For semi-major axis less than two earth radii we exclude all black holes or CDOs with 
masses larger than 6.7 x 10"? kg. Lower mass primordial black holes may be strongly constrained 
by Hawking radiation. We conclude that near Earth black holes are extremely unlikely. 


Primordial black holes (PBH), made by density fluctu- 
ations in the early universe, are an attractive candidate 
for dark matter [I] B]. These PBH could have a large 
range of possible masses depending on the scale of the 
density fluctuations. However, dark matter made of PBH 
with masses Mpgy < 104 kg is constrained by Hawking 
radiation [3] while Mpgy > 10!° kg is constrained by 
microlensing [4}{8]. This leaves a very interesting open 
window of PBH masses between 1014 and 10!° kg that 
is largely unconstrained by present observations. For ex- 
ample, De Luca et al. explain NANOgrav observations 
[9] with gravitational waves from the formation of PBH 
with masses in this range [10]. 


There could be black holes in the solar system. Indeed, 
Scholtz and Unwin speculate that Planet 9, a possible 
outer solar system object of several Earth masses, is a 
PBH [[ [12]. If dark matter is composed of PBH, we 
expect ~ 103(10'4kg/Mppy) PBH to be present in the 
solar system (of volume ~ 10° AU?) at any given time. 
It is possible that one of these objects could approach 
Earth or even be captured into orbit through a three 
body interaction. 

More generally, it is possible that Earth may have one 
or more “dark moon” in addition to Luna. These ob- 
jects could be composed of PBH or other forms of dark 
matter that allows them to escape detection via electro- 
magnetic observations. In addition to PBH, other forms 
of dark matter may concentrate into compact dark ob- 
jects (CDO). These objects are assumed to have small 
non-gravitational interactions with normal matter. Some 
possibilities or names for CDO include Boson Stars [I3], 
Dark Blobs [14], asymmetric dark matter nuggets [I5], 
Exotic Compact Objects [I6], Ultra Compact Mini Ha- 
los (UCMH) [I7] made for example of axions [18], and 
Macros [I9]. 

The lifetime of near Earth orbits for conventional ob- 
jects is short because of air resistance. Dark objects, 
be they PBH or CDO, may have unique long lived near 


Earth orbits because of their small non-gravitational in- 
teractions. In this letter, we directly search for small 
gravitational accelerations (tides) from a PBH or CDO, 
with a 1014 — 101° ke mass, in a near Earth orbit. 

In previous work we searched for gravitational waves 
from CDO merging with neutron stars or orbiting in- 
side the sun [21]. We also searched for small gravitational 
accelerations from CDO moving inside the Earth [22], see 
also [23H25]. In the present letter we significantly extend 
ref. [22] to search for CDO or PBH orbiting partially 
or completely outside the Earth with a range of longer 
orbital periods. 

One might observe electromagnetic radiation as mate- 
rial accretes onto a black hole, see for example [12]. Even 
if there is little accretion, one can still search for a dark 
object’s gravity by accurately tracking spacecraft, see for 
example [27], or using a gravimeter on Earth. 

Gravimeters [28] measure the local acceleration due 
to gravity and can detect small tides from an orbit- 
ing CDO or PBH. Sensitive superconducting gravimeters 
have been deployed at several locations around the world 
[29]. They are used to observe a wide range of geophys- 
ical phenomena including Chandler wobble, solid Earth 
tides, post glacial rebound, seismic free oscillations and 
hydrological processes [30]. In addition to geophysics, 
they have been used to search for a dependence of grav- 
ity on a hypothetical preferred reference frame [31][32}, or 
the violation of Lorentz invariance [33] [34], as the Earth 
translates or rotates. In addition, gravimeters have been 
used to search for oscillations of the Earth excited by 
gravitational waves [35]. 

An object orbiting through or around the Earth with 
coordinate r(t) will have an acceleration, 


r(t) x, Menal?) 
m ~ r. (1) 


r3 


Here G is Newton’s constant and Mene(r) is the enclosed 
mass of the Earth that is interior to r. This reduces to the 


mass of the Earth Mpg for r > Rp the radius of the Earth, 
see ref. for details. We numerically integrate Eq. 
using the Velocity Verlet algorithm with a time step 
At = 1 s, typically for a total time of 27° s. This simple 
procedure works for orbits inside, partially inside, and 
outside the Earth. We assume the Earth is spherically 
symmetric. Explicitly including the Earth’s quadrupole 
moment does not significantly change results for excluded 
masses, see below. We also neglect small perturbations 
from the moon or the sun. These are expected to be 
unimportant for r near the Earth. 

A gravimeter located at R(t) will feel a time-dependent 
acceleration a(t) that is the difference of the gravitational 
acceleration from the compact object minus the acceler- 
ation of Earth towards the object, 


Mene(T) 
r 
Mr r3 


r-R z 
Sl -R. (2 
Here Mp is the mass of the dark object, either Mppy 
or Mcpo, and Ê = R/R. Note that the gravimeter 
only measures the vertical component of a. We choose 
a coordinate system where the gravimeter is located at 
R(t) = Re(cos 0; cos wt, cos 8; sin wt, sin 01). Here 0; is 
the latitude of the gravimeter (48.33° N for the Black 
Forest Observatory) and w = 2r /day. We fast Fourier 
transform (FFT) a(t) from Eq. |2| (including a Hanning 
taper) using the numerically integrated orbit r(t). This 
yields the signal amplitude A( f) for frequency f that will 
be compared to the FFT of gravimeter data. We will be 
interested in f from ~ 107° to 2 x 1074 Hz depending 
on the orbit. 

We now analyze gravimeter data. Data from super- 
conducting gravimeters (SGs), deployed at various lo- 
cations around the world, has been archived by the 
Global Geodynamics Project (GGP, 1997-2015) and by 
the International Geodynamics and Earth Tide Service 
(IGETS, 2015-) [29}[36]. We focus on the SG at the Black 
Forest Observatory (BFO at 48.33°N, 8.33°E) in South 
Western Germany because of its low noise. We also con- 
sider the SG at Djougou, Northern Benin (DJ at 9.74°N, 
1.61°E) in West Africa because of its location near the 
equator. 

We analyze 11.1 years of BFO gravimeter data from 
Dec. 1, 2009 to Jan. 23, 2021 [37] and 7.9 years of DJ 
data from Jan. 1, 2011 to Nov. 30, 2018 [38]. We an- 
alyze this data as in Ref. [22], but now focusing on the 
lower frequency band ~ 107° to 2 x 1074 Hz. We care- 
fully handle artifacts in the data: times when the instru- 
ment behaved non-linearly due to saturation from large 
quakes, operator interference or other malfunctions. We 
subtract from the data a synthetic tidal model for the 
stations that includes the effect of ocean loading [39]. In 
addition we partially correct for accelerations due to at- 
mospheric mass fluctuations above the gravimeter based 
on the barometric pressure recorded near the gravimeter, 
see Ref. [22] for details. 


a(t) = GMp{ 


In Figs. }1] and [2] we show Fourier amplitude spectra 
of the pressure corrected and Hanning tapered gravity 
residuals for BFO and DJ, respectively. These spectra are 
normalized such that a pure time domain signal h(t) = 
Ag sin(27 fot) of amplitude Ag yields A( fo) = Ao. The 
average noise and standard deviation for a number of 
frequency bins is shown along with a conservative upper 
limit that is 100 above the average noise. A signal from 
a PBH above this limit can confidently be ruled out [40]. 
We fit these upper limit noise curves with 


1 ale 3 


Anrolf) = 1.25 x 107 ( 


for the BFO station and 


Apx(f) ~ 4.71 x 10-3(* 7) pm/s’ (4) 


for the DJ station. These fits are valid for 107° Hz < 
f <2 10-* Hz. 
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FIG. 1. Fourier amplitude spectra of gravity residuals for 


the Black Forest Observatory gravimeter versus frequency f. 
Shown in black are the mean and standard deviation ø. The 
red line is a conservative upper limit of mean plus 100 and the 
dashed blue line is the simple fit to this limit given by Eq. 
The three peaks in the upper limit at frequencies above 107 
Hz are solar and lunar diurnal, semidiurnal, and terdiurnal 
residual tides. 


We compare Fourier amplitude spectra A(f) for a 
given orbit to the gravimeter upper limit noise curves 
in Eqs. to deduce excluded masses Mez. The ex- 
cluded mass is the smallest object mass Mp in Eq. [2] so 
that A(f) just equals the noise limit curve for one har- 
monic frequency f. Any object with Mp > Mex can be 
ruled out by the gravimeter data. Table |I| presents ex- 
cluded masses for some selected orbits as shown in Figs. 
The final excluded mass for a given orbit is the 
minimum of the excluded mass based on BFO and DJ 
observations. 
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FIG. 2. Fourier amplitude spectra for the Djougou gravime- 
ter as in Fig. The dashed blue line is the fit to the upper 
limit in Eq. 
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FIG. 3. Fourier amplitude A(f) vs frequency f for equatorial 
(inclination 6; = 0) orbits observed with the DJ gravimeter. 
Eccentricity € = 0 orbits with semimajor axis a = 3Rp (thick 
black curve) and a = 6Rpg (blue curve) are shown in addition 
to an elliptical orbit with a = 6.05Rpx and eccentricity € = 0.1 
(dotted red curve). In each case the object mass is equal to 
the appropriate excluded mass see Table |I| The dashed blue 
curve shows the DJ noise limit from Eq. 


Figure |3| shows A(f) for DJ observations of equato- 
rial orbits (inclination 6; = 0). As the semimajor axis 
a of the orbit increases, A(f) shifts to lower frequencies, 
where the gravimeter noise is larger. This shift is be- 
cause of the larger orbital period. Elliptical orbits (with 
nonzero eccentricity €) have higher harmonics with larger 
amplitudes compared to A(f) for a circular orbit. As a 
result Mexr for € > 0 is expected to be lower than Mex for 
a circular orbit (with the same a). This can be seen in an- 
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FIG. 4. Fourier amplitude A(f) vs frequency f for circular 


orbits of radius a = 3Rg as observed by the BFO gravimeter. 
The thick black lines are for inclination 0; = 0° and the thin 
red lines are for 0r = 45°. In each case the object mass is 
equal to the appropriate excluded mass see Table|I| The blue 
dashed curve shows the BFO noise limit from Eq. 


a (Re) Or (°) e Gravimeter Mex (kg) 


3 0 0 DJ 4.0 x 1014 
6 0 0 DJ 3.7 x 1016 
6.05 0 01 DJ 2.7 x 101° 
3 0 0 BFO  5.0x 10! 
3 45 0 BFO 3.0x 10! 


TABLE I. Excluded mass Mex for example orbits with semi- 
major axis a, inclination 67, and eccentricity € as observed 
with the indicated gravimeter. 


other way. Elliptical orbits come closer to the gravimeter 
than do circular orbits, with the same a. This increases 
the gravity signal of the compact object and leads to a 
smaller Mex. 

Figure [4] shows A(f) for BFO observations of circular 
orbits. As the inclination of the orbit r increases from 
0 to 45°, Mex decreases. In general, inclined orbits will 
come closer to the BFO gravimeter, at latitude 48.33° N, 
than do equatorial orbits. This leads to larger A(f), for 
a given Mp, and a smaller excluded mass Mex. Further- 
more, inclined orbits will have more harmonics in A(f) 
involving sums and differences of the orbital frequency 
and multiples of the Earth’s rotational frequency. This 
is also illustrated in Fig. 

Figure [5] shows excluded mass Mer versus a for equa- 
torial orbits 0r = 0. Here the DJ gravimeter dataset 
provides tighter constraints than the dataset from BFO. 
Furthermore both gravimeters are much more sensitive 
for elliptical orbits than for circular orbits. 

In general Mex increases smoothly with a. However, 
there are resonances where the orbital period over the 
Earth’s rotational period is equal to the ratio of small 
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FIG. 5. Excluded mass Mex versus a for equatorial orbits 
0r = 0. Black circles are for BFO, and red squares DJ, obser- 
vations of circular orbits. Elliptical orbits observed with DJ 
are shown for eccentricities « = 0.25 downward red triangles, 
0.5 diamonds, and 0.75 upward triangles. 
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FIG. 6. Excluded mass Mex versus a for circular orbits. 
The red squares show DJ observations for 6; = 0 while the 
black circles show BFO observations of inclined orbits with 
r > 30°. The dashed black line is Mez ~ 2.1 x 10'?(a/Rz)® 
kg, see Eq. 


integers. Resonances can introduce structure in Meg ver- 
sus a. There is a large peak in Mex for equatorial orbits 
with a ~ 7Rg. This corresponds to nearly geostation- 
ary orbits. Here objects move very slowly with respect 
to the gravimeter. As a result the gravimeter is almost 
blind to their presence. This blind spot is only present 
for nearly equatorial circular orbits and largely vanishes 
by 6; = 15°. 

Excluded masses for a < Rp in Fig. correspond 
to objects moving inside the Earth as discussed in ref. 
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[22]. For these objects there is a second blind spot at the 
center of the Earth a = 0. Here the object is also at rest 
with respect to the gravimeter. 

Figure [6] shows Mex for circular orbits with a range of 
inclinations. An approximate fit to an upper limit for all 
of the Mex in Fig. [6]is 


Mex © 2.1 x 10!2(<_)” kg. (5) 
Reg 
This is shown as a dashed line in Fig. Further- 
more since Mex is smaller for elliptical orbits, Eq. [5] also 
serves as a conservative upper limit for elliptical orbits. 
Therefore Eq. [5] provides an upper limit to Mex that is 
valid for orbits of any eccentricity or inclination. It is 
our primary result. For example, at a = 2Rg we have 
Mex = 6.7 x 1013 kg. A spherical object of this mass and 
terrestrial densities has a radius of a few kilometers while 
a black hole of this mass has a Schwarzschild radius of 
100 fm. 

The a® dependence of Mer in Eq. |5| can be under- 
stood as follows. First, the tidal force falls as ~ 1/a?. 
Second, the orbital frequency f decreases as 1/a?/? (Ke- 
pler’s third law). Finally, the noise limit in Eqs. |3] and 
[4] increases as f decreases. Therefore the ratio of signal 
to noise decreases faster than 1/a? so that the excluded 
mass grows faster than a3. 

Our search is limited by the background in the 
gravimeter signal from atmospheric fluctuations. It 
should be possible to somewhat increase the sensitivity 
with a coherent search involving gravimeters at several 
locations. However, this would likely involve a great in- 
crease in complexity because one would need to analyze 
gravimeter data separately for every compact object or- 
bit. 

Searches around other solar system bodies are also pos- 
sible. The seismometer on the InSight lander on Mars has 
sensitivity as a gravimeter. Unfortunately, this instru- 
ment is exposed to some wind and thermal noise [41]. Ob- 
servations on the Moon could largely avoid atmospheric 
fluctuations. Indeed a gravimeter was deployed during 
Apollo 17 [42] 43] but it did not function well. 

We searched for periodic gravimeter signals from small 
tidal interactions of dark objects and found none. We 
rule out all near Earth objects orbiting with semima- 
jor axis a and mass M > M.,(a), see Eq. We find 
Mex(2Re) © 6.7 x 10'8 kg and rule out all near Earth 
black holes orbiting with a < 2Rg and M > 6.7x10' kg. 
Furthermore, dark matter made of primordial black holes 
of mass Mpgy < 10! kg is significantly constrained by 
the lack of observation of Hawking radiation [3]. There- 
fore a near Earth black hole may be extremely unlikely. 
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